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Abstract.  The  usefulness  of  the  canopy  flow  index  concept  is  demonstrated  for  a  two-story  evergreen  tropica] 
forest.  A  sample  of  about  2500  wind  profiles  was  utilized.  It  encompasses  a  large  range  of  ambient  wind 
conditions  and  spans  the  whole  monsoon  cycle  in  Southeast  Asia. 

It  was  found  that  the  use  of  two  canopy  flow  indices  (one  for  the  upper  and  one  for  the  lower  canopy) 
would  be  necessary  to  simulate  the  average  canopy  flow.  For  the  upper  canopy,  an  average  value  of  4.04  was 
obtained;  for  the  lower  canopy  an  index  of  1.77  was  computed.  The  indices  seem  to  be  independent  of  the 
ambient  wind  speed  (if  2  m  s  '  1  is  exceeded),  yet  strongly  dependent  on  wind  direction. 


1.  Introduction 


The  structure  of  the  turbulent  air  flow  in  vegetative  canopies  is  quite  complex.  The 
turbulence  within  the  stand  is  generated  by  the  breakdown  of  the  atmospheric  flow  into 
eddies  (externally  induced)  and  by  eddies  generated  by  the  many  inhomogeneous 
elements  that  comprise  the  canopy  (internally  induced).  Historically,  the  approach  to  the 
problem  of  canopy  flow  proceeded  along  the  following  directions: 

(a)  theoretical  —  using  basic  equations  of  motion  and  continuity,  supplemented  by 
semi-empirical  parameterizations; 

(b)  semi-empirical  -  based  on  wind  tunnel  studies  and  field  experiments. 

Most  of  the  mathematical  models  assume  steady  state  flow  in  a  horizontally  infinite 
canopy  under  neutral  thermal  stratification.  The  basic  differential  equation  so  derived 
is: 


where 


a, 

dz  2 


K(z)  =  eddy  viscosity, 
v  =  molecular  viscosity, 

U  =  mean  horizontal  wind  speed, 

CD  =  the  form  drag  coefficient, 

A(z)  =  the  leaf  and  branch  area  per  unit  volume 
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The  left  term  represents  the  divergence  of  the  Reynolds  stress  and  the  right  term  is  equal 
to  the  fluid  drag.  This  equation  can  be  solved  by  numerical  techniques  and  the  methods 
used  vary  from  investigator  to  investigator.  Some  use  the  eddy  viscosity  hypothesis  (Tan 
and  Ling,  1963)  to  parameterize  the  transport  of  momentum;  some  prefer  the  mixing- 
length  hypothesis  (Inoue,  1963;  Cowan,  1968;  Cionco,  1965);  some  use  the  second- 
order  closure  technique  (Shaw,  1977)  while  others  use  the  structured  fluids  approach 
(Silbert,  1970).  An  example  of  the  empirical  treatment  is  the  work  of  Uchijama  and 
Wright  (1964)  who  postulated  a  wind  speed  function  which  included  a  term  for  the  effect 
of  the  plant  community  on  wind. 

Considerable  understanding  of  canopy  flow  was  gained  from  controlled  wind-tunnel 
studies  conducted  by  various  investigators  such  as  Plate  and  Quraishi  (1965),  Sadeh 
et  at.  (1971),  Seginer  etal.  (1976).  Wind-tunnel  studies  in  contrast  to  studies  at  natural 
sites  allow  one  to  simplify  real  situations. 

The  present  study  utilizes  the  formulation  of  Cionco  ( 1 965)  who  derived  the  following 
canopy  flow  equation: 

2  01/ 1"  d(ln/c)  +  Pul  _  SIP 

h  dx  lh2dx  dx  +  h23x2\  ~  Pc  (2) 

where 


x  =  z/h , 

z  =  height  above  surface, 

h  =  canopy  height, 

lc  =  mixing  length  in  canopy, 

5  =  1/2  CdA(z). 

U  sing  the  assumptions  of  an  idealized  canopy  with  a  mixing  length  that  is  constant  with 
height,  this  equation  can  be  solved  yielding: 

l/(z)=  (/(/!)  exp  [-a(l- z/h)]  (3) 

where 


a  =  extinction  coefficient, 


(«/A)3  = 


CdA(z) 
21*  ’ 


i.e.,  the  wind  speed  falls  off  exponentially  with  a  depth  scale  proportional  to 
2'li\2l\CDA(Z)Y ,/3,  where  A(z)  and  CD,  or  their  product,  are  assumed  constant  with 
height. 


2.  The  site 

The  experimental  site  is  a  tropical  dry  evergreen  forest  of  about  80  kmJ,  situated  at 
14°31'  N,  101 0 55'  E  in  Thailand  at  elevations  between  300  and  600  m.  The  region  is 
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under  a  monsoonal  regime  with  two  distinct  seasonal  circulations:  a  dry  winter  outflow 
from  a  cold  continental  anticyclone,  and  a  moist  summer  inflow  into  a  continenfal  heat 
low.  The  climate  is  classified  as  tropical  savannah  with  annual  rainfall  of  15o0  mm. 

The  wind  speed  data  utilized  in  this  study  were  collected  along  a  50  m  high  tower.  The 
tower  was  instrumented  for  temperature,  dew  point,  wind  speed  and  direction,  rainfall 
and  radiation  measurements.  Figure  1  illustrates  the  location  of  the  measurement  levels 
of  the  wind  speeds  along  the  tower. 


Fig.  1 .  The  distribution  of  the  wind  speed  measurement  levels  along  the  forest  tower. 


The  wind  speed  sensors  were  Climet  three-cup  anemometers  having  a  threshold 
velocity  of  0.3  m  s~ A  totalizing  signal  was  read  every  10  s.  The  instruments  were 
calibrated  before  and  after  the  experiment  at  the  National  Bureau  of  Standards  (NBS), 
Washington,  D.C.  in  1968  and  1970.  The  uncertainty  of  the  wind  speed  measurements 
prior  to  the  experiment  was  reported  as  follows: 

Air  Speed,  m  s’ 1  >13.4  4.5  2.2  1  0.5 


2  4  6  8  10  meter* 

RADIUS  OF  CROWNS 

Fig.  2.  The  vertical  distribution  of  the  canopy  crown  radii,  after  Cionco  (1980). 

distribution  of  canopy  crown  radii  (Figure  2)  was  constructed  (Cionco,  1 980).  The  tallest 
trees  are  20-35  m  high  with  a  dense  crown  canopy  of  Hopeaferrea;\he  tops  ofttyasecond 

TABLE  1 

Two-week  averaged  wind  velocities  at  selected  levels  along  the  forest  (FT)  and 
clearing  (CT)  towers  (F,  -  at  the  46  m  level,  F2  -  at  the  32  m  level) 


Period 

fin  1Q7A\ 

Level 

CT 

FT 

(in  i  y  iyj) 

m 

Fms'1 

vl/v2 

Fm  s'1 

January 

46 

2.84 

1.17 

2.200 

3.88 

32 

2.42 

0.566 

February 

46 

6.47 

1.18 

4.84 

2.15 

32 

5.47 

2.25 

April 

46 

4.75 

1.17 

3.63 

2.40 

32 

4.06 

1.51 

June 

46 

5.73 

1.15 

4.17 

2.20 

32 

4.98 

1.85 

July 

46 

6.49 

1.14 

4.70 

2.20 

32 

5.67 

Z13 

September 

46 

5.16 

1.08 

3.93 

2.27 

32 

4.77 

1.73 
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story  canopy  are  5— 17  m  tall.  The  species  here  were  Hydrocarpus  iUcifolius,  walsura 
trichostemon,  and  Memecyclon  ovatum.  The  undergrowth  is  composed  of  seedlings  2-3  m 
tall. 


3.  Data 

Wind  speed  data  for  January,  February,  April,  June,  July,  and  September  were  used. 
Approximately  ten  days  of  observations  were  available  for  each  month.  Each  day  had 
48  wind  profiles  averaged  half  hourly;  each  half-hourly  value  was  based  on  180  10-s  data 
points.  The  average  wind  speeds  at  about  20  m  above  the  forest  canopy  are  presented 
in  Table  I.  Similar  information  obtained  along  a  tower  in  a  nearby  clearing  is  also 
included,  to  illustrate  the  drag  effect  of  the  forest  canopy  on  the  ambient  wind  field.  The 
average  wind  direction  frequency  distribution  in  the  region  is  illustrated  in  Figure  3. 
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Fig.  3.  Wind-rote  diagram  for  Nakhon  Ratchasima  situated  60  km  north  from  experimental  site  for 
1956-1960,  from  ASRCT  (1967).  The  number  in  center  represents  percentage  frequency  of  calms. 


During  the  winter  monsoon  (November-February),  the  prevailing  winds  are  from  the 
NE,  while  during  the  summer  monsoons  (May-September)  the  prevailing  winds  are  from 
the  SE;  the  steadiest  SE  flow  occurs  in  June-July. 


4.  Procedures  and  Results 

Before  adopting  a  unified  approach  to  represent  the  canopy  flow,  numerous  profiles  for 
each  season  were  computer  plotted.  It  became  evident  that  the  profiles  fit  into  one  of 
the  following  categories: 

(1)  an  exponential-type  decrease  from  canopy  top  to  the  lowest  measurement  level, 
with  lower  rate  of  decrease  for  z/h  <  0.4  (i.e..  Figure  5a,  July,  1970). 

(2)  an  exponential-type  decrease  down  to  z/h  =  0.6;  an  increase  in  wind  speed  at 
z/h  =  0.4;  a  lower  rate  of  decrease  for  z/h  <  0.4  (i.e..  Figure  5b,  April,  1970). 

a 


Fig.  4a.  Computer  printed  least-square  fit  to  the  computed  canopy  flow  index  a  for  June  1970,  using 
profiles  for  which  the  wind  speed  at  the  highest  measurement  level  exceeded  5ms'1.  The  number  of  cases 
with  equal  a  is  represented  by  a  numeral  X  was  used  when  the  number  of  cases  exceeded  ten. 


ON  THE  -CANOPY  PLOW  INDEX  OF  A  TROPICAL  FOREST 


3 


0  5  C  15  20  25  30  36  40  45  50  2  (m) 

Fig.  4b.  Computer  printed  least-square  fit  to  the  computed  upper  canopy  Row  index  a  for  April  1970.  using 
profiles  for  which  the  wind  speed  at  the  highest  measurement  level  was  between  34  ms  '1. 

(3)  an  exponential-type  decrease  down  to  zjh  =  0.6;  a  no-shear  layer  between 
0.4  5  z/h  i  0.6;  a  lower  rate  of  decrease  for  zjh  <  0.4  (i.e.,  September  1970,  not  illus¬ 
trated). 

These  differences  are  also  evident  from  the  distribution  of  the  canopy  flow  index  a  for 
the  various  cases  (Figures  4a-b).  For  June  (Figure  4a),  the  a  values  cluster  around  two 
‘centers’  (namely,  4.0  and  1.77);  in  April,  due  to  the  frequent  ‘jet’  feature  in  the  canopy 
flow  at  z/h  «  0.4  level,  the  a’s  are  less  uniform  in  magnitude,  resulting  in  high  standard 
deviations.  Consequently,  the  following  procedure  was  adopted. 

For  each  period  of  study,  the  data  were  stratified  according  to  wind  speed  at  the  highest 
measurement  level  above  the  canopy;  the  attenuation  coefficient  a  was  computed  for 
each  such  group  separately  using  Equation  (3).  A  least-square  fit  was  obtained  for  the 
computed  a  values.  Table  11  presents  summaries  of  alfthe  results,  i.e.,  average  a  for  upper 


TABLE  11 

Statistical  summary  about  the  computed  canopy  flow  index  a 


O  B 

1  T 

T3 


in  O  ‘n  N  on  n  -  o  O'  o  di  oo  o  N  n  no  >6 

I  I  i  l  l  ~  ~  —  do  -Joood  o  —  d  d  d  —  —  —  —  d  o 


g  g  - *  rw  O  fN  N  vO  >C  O  mmromg  ^  °  ^  o  ^ 

I  i  i  |  i  ddddo  ddddd  ddodd  odd  od  o 


r-~  r-'  oo  oo  r-  «  «  «  in  ^  oo  oo  oo  qo  n  oooooooooo  ofe 

I  f  I  I  l  ddodd  ddddd  ddddd  ddddd  d 


O  t  o  n  -  Q  O  O  r^i  «  n  -  Ti  O'  OO  o  <N  —  o 
^  ^  rw  <n  <s  ^  ao  so  ncnm<Nr^  n  w  »N  n  oo  rn 
I  I  l  I  l  ddddd  ddddd  ddddd  ddddd  d 


m  r*  r>  ©  r~  ir>  —  —  r-  —  Q  O  oo  m  m  «/■%  sn  sn  in 

-o  «"»  ■*»•  «/■>  o  m  m  n  m  <s  miaow  r^r-oonoo  m 

I  l  I  I  I  ddddo  -  -  n  wwi  d  d  d  o  ~  ddddd  ~ 


O  r-  os  —  o^r^m  m  <n  r«  n  is  m  o  »  o  so 
OOiwN  x  35  i  »n  -  ic  Xh  «  ^  p-oooooor-  d 

i  i  i  i  i  H  ri  -  -  ^  ~  ^  ^  ;  ~  ~  ;  —  — ;  —  —  — ;  — 


8  0  o«  m  t  ro  *r>  —  «r>  —  «o  —  rw  >o  r-  oo  ©  «—  O' 

—  »  O  «  «  W  N  If  h  M  n  N  o  W  fS  f»0  so  00  O  vO 

I  I  I  I  I  Tf  rn  ^  rn  ^  ^  w  H  ^  ^  sf  n  r!  i*i  n  w  ri  <n 


c  6  w 

5  8  s 


s||  | g |  ? 00  assas  ^ a  ^ -  si  2  §  5  5  s 


4.16  1.62  1.99  0.40  0.77  0.16  0.26 

4.45  1.67  2.20  0.35  0.78  0.18  0.12 

4.97  1.60  2.82  0.32  0.78  0.22  -  0.31 

4.27  1.39  2.31  0.33  0.79  0.18  -0.29 
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and  lower  canopy;  correlation  coefficients;  intercepts  and  slopes.  As  evident  from 
Table  II,  the  highest  frequency  of  wind  speeds  in  January  is  in  the  low  range  of 
0-2.0  ms1.  Thompson  and  Pinker  ( 1975)  analyzed  the  stability  conditions  for  this  site 
and  showed  that  January  was  the  most  unstable  period  and  as  such,  least  suitable  to 
the  application  of  a  canopy  model  based  on  neutral  stratification.  Therefore,  no  as  are 
presented  for  this  period.  For  all  the  months  treated,  there  does  not  seem  to  be  a 
dependence  of  a  on  ambient  wind  speed  as  long  as  winds  of  2  m  s' 1  are  exceeded.  The 
average  value  of  a  for  the  lower  canopy  seems  to  be  half  of  that  for  upper  canopy.  The 
standard  deviations  seem  to  be  related  to  the  variability  of  wind  direction.  As  illustrated 
in  Figure  3,  April  and  September  are  among  the  most  variable  in  wind  direction  and  the 
standard  deviations  in  a  during  these  periods  are  3  times  those  for  the  other  months. 
Strong  dependence  on  wind  direction  was  also  reported  by  Cionco  (1978). 

The  correlation  coefficients  of  the  regression  equations  are  above  0.8  except  for  April 
and  for  the  low  wind  speed  range  during  other  months.  The  average  values  for  all 
periods  studied  are  summarized  in  Table  III.  The  range  of  the  standard  deviations 
during  the  various  periods  is  indicated.  The  values  of  a  as  compiled  by  Cionco  (1978) 
are  summarized  for  comparison. 

TABLE  III 

Average  canopy  flow  index  a  for  forest  canopies  after  Cionco  (1978)  and  Shinn 
(1971)  and  for  present  study 


Forest  canopy 

a 

mean  ±  st.  dev. 

Gum-maple 

4.42  ±  1.05 

Maple-fir 

4.03  ±0.65 

Jungle 

3.84  ±  1.52 

Spruce 

2.74  ±  1.29 

Oak-gum 

2.68  ±0.66 

Tropical-evergreen 

mean  ±  st.  dev.  range 

upper  canopy 

lower  canopy 

►  5.0* 

3.98  ±  (0.25-1.32) 

177  ±  (0.25-0.40) 

4.0-5.0 

4.05  ±(0.52-1.95) 

1.81  ±(0.30-0.44) 

3.0-4.0 

4.11  ±(0.47-3.51) 

1.72  ±  (0.30-0.70) 

2.0-3.0 

4.06  ±(0.45-3.51) 

1.57  ±  (0.23-0.80) 

0.O-2.0 

3.4  ±(0.67-3.25) 

1.40  ±(0.3 1-0.60) 

*  Wind  speed  range  20  m  above  canopy. 


To  test  the  suitability  of  the  computed  canopy  flow  index  a  to  actual  conditions, 
average  profiles  for  each  period  investigated  were  compared  to  the  model  profile  using 
the  appropriate  canopy  index.  The  results  are  illustrated  for  selected  cases  in  Figures  5 
and  6. 

Figures  5a-b  illustrate  a  normalized  daily  average  profile  on  a  log-log  scale  for  June 
and  April  and  the  corresponding  predicted  profiles  using  the  appropriate  a  value  from 
Table  II.  Figure  6  represents  an  all-season  average  and  the  scatter  obtained  with  a  ±  oa. 
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CANOPY  WIND  PROFILES 
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Fig.  6.  \  seasonal  daily  average  normalized  wind  profile  on  a  log-log  scale  (v).  the  predicted  profile  (A), 
and  the  scatter  around  the  predicted  mean  obtained  with  au,  a,,  and  a„  ±  <t,  +  o,  (V).  (□).  as  given 

in  Table  II. 


The  largest  discrepancies  were  found  in  the  upper  canopy  in  April  and  September  (same 
level)  and  might  be  related  to  the  large  wind  direction  variability  during  these  two  months. 

5.  Summary 

The  canopy  flow  in  a  two-story  evergreen  tropical  forest  was  investigated  using  a  large 
sample  of  wind  profiles  which  span  the  whole  monsoon  cycle  in  Southeast  Asia.  Each 
of  the  wind  profiles  was  an  half-hourly  average  based  on  180  data  values  for  each 
measurement  level.  Inspection  of  the  profiles  led  to  the  adoption  of  the  canopy  flow  index 
concept  as  a  possible  way  to  model  the  canopy  flow.  Numerous  experiments  were 
conducted  to  identify  the  most  appropriate  model  for  this  two-story  canopy.  It  was  found 
suitable  to  use  two  indices;  an  average  value  of  4.04  was  found  for  the  upper  canopy 
(0.6  <z/h<  1)  and  an  average  value  of  1.77  for  the  lower  canopy  (0  <  z/h  <  0.4).  The 
comparison  of  the  modeled  profiles  with  the  measured  averaged  profiles,  seemed  in  good 
agreement.  The  computed  a  for  the  upper  canopy  is  close  to  the  canopy  index  reported 
for  Gum-maple  (4.42  ±  1.05)  and  Maple-fir  (4.03  ±  0.65)  forests. 
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Due  to  the  large  number  of  profiles  available,  this  study  could  address  one  of  the  still 
unresolved  issues  regarding  canopy  flow,  namely,  the  dependence  of  a  on  the  ambient 
wind  speed.  No  dependence  seemed  evident  if  winds  of  2  m  s  1  above  the  canopy  were 
exceeded.  However,  dependence  on  wind  directions  seems  to  be  a  plausible  reason  for 
the  large  variability  in  the  index  as  indicated  by  the  increased  standard  deviation  of  a 
for  periods  of  less  steady  ambient  flow. 
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